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Abstract 

We calculate up to four loops the free energy of the two-dimensional (2D) 0{n) nonlinear cr- 
model regularized on the lattice with the 0-loop and 1-loop Symanzik improved actions. An effective 
coupling constant based on this calculation is defined. 
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Numerical simulations on the lattice are a convenient method to compute expectation values of dimensionful oper- 
ators in a field theory. To convert the expectation value obtained from the simulation into a matrix element written 
in physical units we must fix the lattice spacing a. In the present paper we will focus our attention on the case of 
psj ' asymptotically free field theories. 

^ [ The lattice scale a is usually fixed through the integration of the beta function of the theory. This function can be 

C^l . evaluated in perturbation theory in which case a is written also as a weak coupling expansion (starting from three 
loops) . It is said that the asymptotic scaling regime has been attained if the first universal terms in the perturbative 
expansion of the beta function are enough to have a good knowledge of the lattice spacing a. In principle we may 
expect that the better the weak expansion of the beta function is known, the more accurate the measurement of a 
Qs [ dimensionful physical observable will be. However this hope usually fails due to the bad convergence properties of 
^\ ' the weak expansion. In particular, the asymptotic scaling regime is barely achieved. 

A cure to this problem may come from a redefinition of the expansion parameter in the perturbative series [2J-§| . 

Ci . Usually the expansion parameter is the bare coupling constant g. However we may define any other parameter gE 

I ■ related to the old one by an expansion gE = g + 0{g^). If gs can be numerically determined by some non-perturbative 

Mh, procedure then the beta function, expressed as a power series in terms of gE, can be regarded as a resummation where 

^ ' the non-perturbative effects have been absorbed in the very definition of the new parameter gE- The problem is the 

choice of such a non-perturbative procedure that would produce a good resummation and a rapid convergence of the 

power series in gE- The calculation scheme where g is substituted by gE shall hereafter be called "effective scheme". 

In Ref. W the 2D 0{n) nonlinear cr-model has been analysed. This is an asymptotically free field theory, which 

shares several physical properties with four-dimensional (4D) Yang-Mills theory, among others the spontaneous 

C^ ' generation of a mass in addition to the asymptotic freedom itself. In 0] the authors argued that the density of action 

(the internal energy) is a good operator to define gE (actually the choice of this operator was proposed in the seminal 

papers ||l|,|| for 4D gauge theories as well as for the 2D nonlinear cr-model). In particular the corrections to asymptotic 

scaling in this effective scheme vanish in the large n limit. In Ref. 0, also devoted to the study of the 2D nonlinear 

cr-model, two different local operators were used to define the new expansion parameter and in both cases a clear 

improvement was observed. Moreover the results obtained for the measured quantities in the two effective schemes 

were compatible with each other, which means that the implicit resummations of the series were equivalent. In the 

case of gauge theories other types of redefinitions of the expansion parameter can be used, see for instance [p|. 

Nevertheless in many cases the lack of asymptotic scaling does not suffice to explain the poor results obtained from a 
Monte Carlo simulation. The lattice regularization of any operator (for example the action) differs from its continuum 
counterpart by terms of higher order in the lattice spacing a. If these terms are sizeable then the scaling ratios of 
operator matrix elements as evaluated on the lattice do not behave as prescribed by the continuum Renormalization 
Group equations, i.e.: they do not show physical scaling. This is another source of systematic errors in the Monte 
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Carlo determination of any dimensionful observable. This problem can be mitigated if the Monte Carlo simulation is 
performed with a Symanzik improved action. 

This difficulty occurs again in the 2D nonlinear cr-model. For instance, the calculation of the mass gap in this 
model with 0(3) symmetry has yielded results |7|-[lq] in conflict with the exact analytical value ruHMN of Ref. [[l6[ . 
In Ref. |l^ an energy-based effective scheme was used to extract the mass gap but the deviation between the Monte 
Carlo and exact results was still beyond f 0%. A definite improvement was obtained only after the problem concerning 
the lack of physical scaling was also tackled. In Ref. [g a combination of an effective scheme together with a Symanzik 
improved action allowed an agreement between the Monte Carlo and the exact results well within 2-3% 

Therefore one of the main conclusions of Ref. [p| was that the combination of an effective scheme together with a 
Symanzik improved action may allow to obtain results reasonably clean from all lattice artifacts (assuming that the 
lattice volume is large enough). Our purpose is to further check this conclusion. To this end we want to calculate 
the mass gap of the 2D 0(3) nonlinear cr-model on the lattice by simulating two improved actions, the tree-order 
Symanzik and the 1-loop Symanzik actions ||l^,|l9|. Let us call Am the difference between the mass extracted from 
a Monte Carlo simulation and the exact result of Ref. [^ . We want to study the dependence of Am on the level of 
improvement in the action as well as on the number of corrective terms to the asymptotic scaling in the perturbative 
expansion of the lattice spacing. Firstly we need to calculate these corrections in perturbation theory and this is done 
in Ref. [Q where we calculate the lattice beta function /3-^ for the two above-mentioned Symanzik actions up to four 
loops in terms of the bare coupling g. In the present paper we will show the result for f3^ in terms of an effective 
coupling gE defined through the internal energy E. We will perform the calculation for a generic symmetry group 
0(n), the specialization to the case n = 3 being trivial. 

Our interest is actually twofold because a high precision determination of the mass in the model would also contribute 
to settle a long-standing debate about the validity of the tuhmn value |2|] . 

In the following section we give the expressions of the two Symanzik improved actions and the explicit form of the 
operators for E. We also briefly describe our calculation procedure. In the third section we introduce some notation 
and a list of useful identities that we have used. The final results for the internal energy E will be given in section 
4 together with the calculation of the lattice beta function in terms oi gs- A summary of the results and a brief 
account of the checks performed on them is given in the conclusions. Several appendices are devoted to the finite 
integrals that appear in the final expressions; we list these integrals and their numerical values, explaining how they 
were evaluated and showing some identities which relate them. 

II. THE CALCULATION 

In this section we give an outline of the calculation method that we have followed. The action of the nonlinear 
0(n)-symmetric cr-model on the lattice can be written in many different ways. All of them share the same naive 
continuum limit, i.e. the action of the model in the continuum 

^continuum = 1 I d\ ^ (djix))' , (2.1) 

where g is the coupling constant and is a scalar field with n components constrained by the condition 
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= 1 for aU X . (2.2) 



In this work we shall consider three actions on the lattice. The standard action is the simplest one but yields the 
poorest results in simulations; on the other hand the Symanzik improved actions allow a progressive elimination of 
lattice artifacts as powers of the lattice spacing a and its logarithm log a |19,M . These actions are 



S^'^'^"^^-^ = -J2lk^) ■ Ki ■ ${x) , 
a ^ — ^ 2 

^O-Symanzik ^ ^ ^ |^2^(^) . ^^ . ^^^^ _ _^^(^) . ^^ . ^(^A 
" X ^ ' 

^1-Symanzik = ^ ^ [ i0(a;) • K^ ■ ^(x) - o" C, g (i^l • ${x)) ^ - «' (^6 3 - ^) ^ {^t ^ , 






(2.3) 

where the superscript i-Symanzik denotes the i-loop Symanzik improved action. The several operators in this equation 
are defined by 

Ki ■ ${x) ^^J2 {^^^^^ - "^(^ + ^) " '^(^ - ^^) ' 
K2 ■ ${x) = —Y^ (2${x) ~${x + 2fl)- 4>{x - 2/1) 
1 



9+0(x) 



(x + /2) - 0(x) 



d„ (hix 



a V 



(j)[x - /2) 



(2.4) 



The set of coefficients {ci}i=5^...^9 is determined by the Symanzik improvement program at one loop ||19|,P2| (there is a 
discrepancy between the numerical values of cg reported in Q and ||l^, we agree with the latter, see |23|]). We shall 
often refer to these actions as standard action, 0-loop action and 1-loop action respectively. 

The 1-loop action can be written in a more convenient form (mainly for the Monte Carlo simulation) after intro- 
ducing the operators of Eq. (|2!^) into Eq. (p^) , 
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2 ^ ${x) ■ ${x + 2/2) - 16 ^ (fix) ■ ${x + At) + 

2 ^ ^(x + /2) • (^(x + P) + 2 ^ ^(x) • 0(a: + /2 + P) 

2 ^ 0(a;) • ${x + 2/2) - 8 ^ (^(x) ■ (f{x + /2) 

- 16 ^ ${x) ■${x + Jl) + 2Y, {$ix) ■ ${x + ii)Y + 
2 ^ (0(x) • ^(x + /2)) ((^(x + /2) ■ ${x + 2/2)) + 
Y^ ( (^(a;) ■ 0(x + /2)) (0(x) • ^(x + 9)) + (^(x) • ${x + /2)) (^(x) • ${x - P)) 

(^(^(x) • ${x - Jl)j [${x) ■ 0(x + P) j + [${x) ■ ${x - p)j (^0(x) • 0(x - 17) 

^ (0(x) • <^(x + /2)) - 8 ^ ${x) • 0(x + /2) + 2 ^ (0(x) • 0(x + /2)) (^(x) • ${x - /2) 
M MP 

J2 {Hx + fi) ■ ${x - /2)) - 8^0(x) • 0(x + 2/2) + 
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Y^ I ((^(x + /2) • ${x + P) j + [${x + /2) • ,?(x - P) j + 
((/>(x — /2) • 0(x + i2)j + ((/)(x — /2) • 0(x - i2)j 1 + 



( 4>{x + fi) ■ 4>{x - i^) ] ( 0(a; - p-) ■ (t){x + v) 

( (/)(x + /2) • 4}{x + P) j ( (j){x + ^) ■ (l){x - v) 

{(t>{x + fi)-(t>{x + v)\ (0(a; - /I) • (?!)(x + 9) 

[(j){x - /2) • (p{x - P)j \4)(x + /2) • ^(x - 9) j - 



(2.5) 



We notice that this ex pres sion differs from the one shown in Table 2 of Ref. p4| . We can only say that we have 



carefully checked our Eq.(2^) and that it coincides with the action given for example in Ref. [ pSj . 

The terms proportional to the coefficients Ci are of order 0[g). Therefore the 1-loop action is equal to the 0-loop 
action plus a sum of terms of higher order in g. This fact allows us to compute any perturbative quantity at fc loops 
for the 1-loop action as the sum of the analogous quantity for the 0-loop action plus a set of diagrams with at most 
(fc — 1) loops. 



To the actions showed in Eqs.( |2.3| ) and ( p.5| ) we must still add another term. The constraint shown in Eq.(2.2) is 
introduced under the path integral representation of the theory as a Dirac delta. Hence the partition function Z is 
written as 

Z ^ I Vl{x) n 6 {U(x)) ' - 1) e-^ , (2.6) 

where S indicates any of the lattice actions that we consider in the present paper. This constraint can be solved by 
rewriting the n-component field (^ in terms of a new (n — l)-component field 7? in the following way 

^{x) = {4>i{x), ..., 4>n{x)] — > (7ri(a;), ..., 7r„_i(x), v^l - 7i'i(a;)2 - ... 7r„_i(a;)2 j 

^ {nix), V^-A^r) • (2.7) 

This transformation introduces a change of variables in the functional integration whose jacobian becomes a new term 
to be added to the original action as a measure action |2g] , 

5 ^ 5 + ^ log Vl-7r(x)2 . (2.8) 

X 

Perturbation theory is developed for the (n — l)-component field 7f(a;) around the trivial configuration 7r(a;) = 
for all X. The contribution to perturbative expansions from the measure term vanishes in the continuum when the 
divergences are dimensionally regularized p7| . 

In this paper we will compute the average density of action (or internal energy) E, up to fourth order in perturbation 
theory for all three actions on the lattice. We will begin by computing the free energy F to four loops, from which 
the internal energy reads 

E = ^o + l^f, F.^hm^i^logZ, (2.9) 

where L^ is the volume of a two-dimensional square lattice of side length L and Z is the partition function. The 
constant wo is determined in such a way that the operators representing E are (there is no summation in the index 

^standard ^ ^^^q) . ^(q _^ ^)^ ^ 
^0-Symanzik ^ ^4 .^^^ .^^ ^ ^^ _ 1_^^^^ .^^ ^ ^^^^ ^ 



Then ^ygtandard _ -^ £qj. ^j^g standard action and Wq~ y™™^' = yj^~ y™anzi _ ^^^^^2 for the two Symanzik actions. 

In perturbation theory F is computed through the sum of the connected Feynman diagrams shown in Figs. 1-6. 
The analytic expressions of these diagrams are plagued with infrared (IR)-divergent integrals. An usual method to 
deal with these IR divergences is the introduction of an external source h coupled to the scalar field ]2^,^. We 
have chosen a different IR regulator: we have put the model into a square box of size L with periodic conditions on 
the boundaries. Then we have applied finite-size perturbation theory following Ref. |2|]. This procedure has two 
consequences: firstly all zero modes can be excluded from the sums over momenta and secondly the action contains 
a new term which comes from a Faddeev-Popov determinant, 



(y^v^-A^A 



Sfp = - (n - 1) log I ^ Vl - A^)^ 1 ■ (2.11) 

We have explicitly checked that the diagrams containing vertices from Spp give rise to contributions which vanish 
in the thermodinamic limit L —> oo order by order up to four loops (individual diagrams give rise to non-vanishing 
terms but they cancel when adding up all diagrams at each order). We have not included this list of diagrams in 
Figs. 1-6 because in this paper we give the results for infinite-size lattices (in any case their contribution is rather 
small: for L = 100 it amounts to 0(10"^) ). 

The diagrams in Figs. 1-6 lead to finite sums over internal momenta. Some of the sums become IR-divergent 
integrals in the limit L —> oo. We have worked out these sums in order to separate their finite contribution from the 
IR-divergent one. In order to isolate these divergences we have performed only algebraic manipulations. We have 
avoided other kinds of manipulations, like applying derivatives to the integrands. Then the divergent pieces exactly 
cancel among themselves order by order. At the end, only the convergent sums remain (sums which in the limit 
L —f oo turn into IR-finite integrals) and we can safely remove the IR regulator by sending L to infinity. These finite 
integrals are then calculated as explained in Appendix C. 

The calculation up to four loops for the standard action is known in the literature |g^,^,^,Q. We have checked it 
obtaining analytical results in agreement. The calculation for the 0-loop Symanzik action is known only up to third 
order |3^,|§1 and we have also successfully checked it. Therefore the new results of the present paper are the fourth 
order for the 0-loop Symanzik action and the full calculation to all orders, up to four loops, for the 1-loop improved 
Symanzik action. For completeness, in section 4 we will give the final results for all actions. 

III. NOTATION AND IDENTITIES 

In this section we will introduce some standard notation for the lattice perturbative calculations. Unless otherwise 
stated, we set the lattice spacing a = 1. The sine of half of the /i component of a momentum p is denoted by 

p^ = 2sin^. (3.1) 

The inverse propagator in the standard action is jp = J^nPu- ^'^^ ^^'^ Symanzik action we introduce the notation 

Qp^E^t' (3-2) 

and then the inverse Symanzik propagator, denoted by lip, reads 

^p^F' + ^□p^E^^ 

^'p^Pl + ^Pt- (3-3) 

This propagator is the same for the two Symanzik actions. 

Several identities among momenta are helpful to separate the finite and divergent contributions to E. One of these 
identities involve standard propagators 0, 

_^^_ 2 2 2 ^^ ^^ 

{p + q) + {p + k) + (p + r) =p^ + (f + k"^ + r^ - Spgfer , ^pqkr = E PijMp.kpXp , (3.4) 



and it is valid ifp + q + k + r = 0. For the calculation using the Symanzik actions other relationships are needed, for 
instance [^, 
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p+q 



+ Up+k + n. 



.p^r 



■^pqkr 






11 11 



^sinp^ sing^ sinfc^ sinr^ j 
p / 



(3.5) 



which again requires that p + q + k + r = Q. Moreover, in the calculation of tadpole diagrams we have used 



p + q =p +q- 

2 



~^ —p^cf + odd terms, 
4 

1 



{pTq') = iff + iff + i [if]' {ff - [{ff - Dp) □, - {{ff □,) □ 



-2 P 



-□, 



9-7°. 



n. 



p+q 



n„ + n, -— npD,-— n,Dp 



ip 



3 <? - 2 ^^ '' * ~ 2 ^^ '' ^ ^ ° terms, 



— — npDg + odd terms, 



D 



P+q 



Dp + D^ - -U-pUq - -HgOp + SUpUg + T^^p^q + odd terms. 



(3.6) 



The Feynman diagrams necessary for the calculation are shown in Figures 1-6. 

Let us show two partial calculations as examples of the procedure we have followed. On finite lattices of side length 
L any component of a momentum p can take L discrete values, for instance the first component pi — 2Trii/L {ii = 0, 
1, ..., L — 1). Therefore the sums over momenta are 



L-l 



i-i 



IS- IE-. 

^1=0 ^2=0 



and become integrals in the limit L —>■ oo, 
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d^p 

rr (27r)^ 



(3.7) 



(3.8) 



In Eq. (|3.7| ) the zero mode ^i = ^2 = must be excluded as prescribed in Ref. |2^ . This is the meaning of the stars 
in Eq. (|3.7| ). In the following we use the shorthand 



Z^n ' 



(3.9) 



to denote the sum in Eq.(3.7) when it sums over the momentum p. Besides, although at finite L we are dealing with 
discrete sums, we will often call them "integrals" and the expression summed will often be named "integrand" . 

Momentum conservation is expressed through a Kronecker delta, L'^ S'^{p + q + k + ■ ■ ■). p, q, fc, ... are momenta 
which satisfy pi — 2tt£jJL, etc. The argument of the delta function is assumed to be periodic modulus 2tt. The sum 
in Eq.(3.£), or in Eq.(^), acting on a delta function leads to 



= l--r L^S\q + k + ---) . 



Y L^S^{p + q + k + 
The additional 0(1/L^) contribution in the r.h.s. of this equation can produce terms like 



-T ^- 



(3.10) 



(3.11) 



which are finite after the removal of the IR-regulator. These terms are exclusive of the finite size L regularization: 
notice for example that they cannot be expressed as usual integrals, neither after the L ^ 00 limit. We have checked 
that such terms cancel out when we sum up the contributions from all diagrams at each or der. 

Terms like those in Eq.(3.11), and similar ones coming from the Faddeev-Popov action (2.11), yield finite contri- 
butions to the final result in the ID 0{n) model [E9l. In the 2D model such terms cancel out at least up to fourth 
order. 



To the diagram 6 of Fig ^ it contributes the coefRcient cg of the 1-loop improved action through the vertex 

r» E, E, E. E, E, E. ^^ ^^(^ +'+'=+'-+•'+ ') 

X ^^7r''(p)7r''(g)7r''(fc)7r''(r)7r"(s)7r"(t) sinp,, sing^ sin(fc + r)^ sin(s + i)^ , (3.12) 

Aif abc 

where p, q, k, r, s and t are the momenta which satisfy pi = 2T:ii/ L, etc. and a, b and c the 0{n) indices (running 
from 1 to n — 1). To simphfy the notation we use the same symbol to indicate the field in coordinate space Tf{x) and 
in momentum space 7f(p). The contraction of the six legs to produce three tadpoles can be done in three topologically 
non equivalent ways and leads to 
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E <"-i)"E.,E„E 



sinp^ sinpi, sin {q + k) sin (q + k) ^ 



^v 



^ ' ^-p ^q ^k 

2("-i)E.E,E 



V ^-^1 ^-^k XlpIIqllfe 

sinp^ sin q^ sin {q + k)^ sin {p + k)^ 
sinpp sin q^, sin {p + k) sin {q + k)^ 



p ^ — ^q ^ — ^k 



UpHgUk 



(3.13) 



The firs t integral is non zero only ii fi — v and, by using sin p^ — pf^ — 1/4 pt and the first and last identities in 
Eq.(^, gives 
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eg g^ - {n-iy 1 - -Yi 2Zi - 1 + -Yi + O 
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(3.14) 



where Yi in the limit L ^ oo are finite integrals defined and calculated in Appendix A and Zi are the above-mentioned 
sums which in the thermodynamic limit diverge and that at the end of the calculation must disappear, 



^p { n, 



(3.15) 

Notice that this expression is well defined as long as L is finite because the zero mode lip = is missing. It becomes 

an ill-defined integral only in th e the rmodynamic limit. For large L we have that Zi ~ logL, Z2 ~ L^, Z3 ~ L*, etc. 

The second integration in Eq.(3.13) yields, after some algebra and taking into account again that it vanishes when 



3 / IN Y- V- V- E^sin^PM sin^^P sin^ ^^ 



^p ^-^^q ' ^k 



Hp n, Hfc 



eg 5^ - (71 - 1) 1 - -Fi 2Zi - 1 + -Fi + O 



logL 



(3.16) 



Finally the third integral in Eq.( 3.13 ) needs some algebra to eliminate several odd terms like EpSi'^-PA'/np- After 
this work, it can be rewritten as 



2cg/(n-l)VE. E 



p ^ <q 



X^^j^sin^P/x sin^g^ cosfc^cosfci. 

Hp n, Hfe 



i,,\V.^ 5 1, 



Cg g-" - in-l) 1 - -Yi 4Zi - - + -Yi + Y2 1 

9^2^ ' \ 3 / V 4 6 576 ' 
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(3.17) 



where the property cosA:^ = 1 — 1/2 fc^ was used. 

The second example of calculation is taken from the coefficient of cg in the diagram 5 of Fig. |6[ The vertices with 
four legs are 



i^'^»E,E,E.E.^^*^(' +'+'•+'-) 

X yTT''{p)Tr''{q)TT\k)TT\r 



p + k 



p - k ) (q - r )-2 {p + k ] (q - r 



(3.18) 



from the black spot and 



X j;i''(p)i"(,)i'(t)i'(,-)n,+, , 



ab 



which is the vertex coming from the 0-loop Symanzik part of the action. The contraction gives 



g — Co 



{n-iyh + in-l){h+l2) 



(3.19) 



(3.20) 



where 



^i-E„E„E.E.^^^^(^+'^ + ^+^) 



n 



p+q 



(pTfc') + (p'^fc') {o^r^) - 2 (p^k^) {l^r^) 



Up n, Hfc n. 



^^-E„E„E.E.^^^^(^+'^ + ^+-) 



n 



p+q 



(pTg") + (p^q^) [k^r^) - 2 (pTq^) (fc^r') 



Up n, Hfc n. 



(3.21) 



The numerators can be easily worked out. For instance, the square brackets in the numerator of /i can be rewritten 
by using the Kronecker delta in Eq.( ^.21 ) (Ap,g is introduced in Appendix A), 



p+^k) + (p^k ) icf^r^^ - 2 (f+k ) (cf^r^ 



+ (Ap _fc + p2 + fc2^ (a,^_^ + q^ +?' 



^2 I ^2 



-2 (Ap,fc+p2 + A?2) (Aq,_^ 

— ^p,k ^q,r r L\p^ — k ^q,-~r ~^ ^ ^p,k ^q,~r 



(3.22) 



which is true under the integration. After taking the limit L — > cx), this expression leads immediately to the final 
result for /i — Su + S14 — 2 S15. Analogously I2 = S'lo + •S'le — 2 S'17. The integrals Si are defined and evaluated in 
Appendix A. 
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FIG. 1. Feynman diagrams contributing to the free energy of the standard and 0-loop Symanzik actions at two loops. The 

lines represent the propagation of the scalar field tt, the cross stands for a vertex from the measure action. 
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2. The same as Fig. hi at three loops. 
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FIG. 3. The same as Fig. hi at four loops. 
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FIG. 4. The Feynman diagrams necessary for the calculation of the free energy at two loops for the 1-loop Symanzik action 
are those shown in Fig. hi plus the diagram displayed in this figure. The black spot indicates a vertex proportional to some 
coefficient Ci. 
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FIG. 5. The free energy at three loops for the 1-loop Symanzik action is calculated by adding up the diagrams of Fig. H to 
those displayed in this figure. Same notation as in Fig. H. 
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FIG. 6. The free energy at four loops for the 1-loop Symanzik action is calculated by adding up the diagrams of Fig. B^ to 
those displayed in this figure. Same notation as in Fig. H. 



IV. RESULTS FOR E AND EFFECTIVE SCHEME 



In this section we will give the results for the internal energy E and will explain how to obtain the coefficients for 
P^ in the corresponding effective scheme. We will show the results for all three actions. Those for the standard action 
are not new [|o|, |j3l|g but we have checked all of them. We write the perturbative expansion of E as 



E = Wo — wi g — W2 g 



W3 g^ - Wi g'^ 



(4.1) 
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The integrals -K, J and -ff^ for i = 1, . . . ,5 are defined and calculated in ||4|, p[. They are shown in our Appendix A. 

The results for the 0-loop Symanzik action are known up to three loops [p2lp| and here we have added the fourth 
order. The first terms Wq~ y™^"^' ^ ^ yj^~ ymanzi j-^^^^g been checked obtaining full agreement. The complete set of 
coefficients is 
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O-Symanzik ^ („ _ i) [ _ A _ l^g _ Ai/f _ li/S _ i^S ^ ]_jS _ i JS + lj5 

^ ^ M 256 96 16 ^ 32 2 8 ^ 16 8 16 



S3 y3 91X15 , 5 91 535 , 
4096 5308416 1530 18432 221184 

Yi Y2 Y2 -\ Yi Y2 Yi Y2 

2654208 55296 331776 7962624 

^ -Ys ^Y, Y, + -^Y,' Y, - -^ri^ ^3^ 



6912 27648 331776 11943936 

+ („ _ 1)2 (±_ + 1^ + Lh^ + LhI + -Hi + -K^ --KS + A^ 
^ M 512 96 16 ^ 32 2 8 ^ 16 8 16 



in -if 
32 



32 ^8 ^ 1536 1152 ^ ^ 16 32 ^ 

±Tsy^ + JA^^. _ ^^^^3 ^ 2725 \ 

64 12288 221184 10616832 j 

Hi . (4.3) 



The integrals K^ , J^ , K^, J^ , K^, J^ , Yi, Y2 and Yz were introduced in g. They are shown for completeness 
in Appendix A. The new integrals are 6*9, ijf , iJf , i?^, Hf and H^ . All of them are listed and evaluated in the 
Appendix. 

The expression of the internal energy for the 1-loop Symanzik improved action is a new result of the present paper 
for all loops and is written in terms of the following coefficients 

1 — Symanzik — Symanzik 
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w 



1 — Syinanzik 



— Symanzik 



1 — Syinanzik 



W 



1 — Syinanzik 



— Symanzik 



— Symanzik 



+ (n-l 



2 

+ {n-lf 

3 



Cb 3 



144^ 



Y2.1 +ceYi 



ci + -C8 



7>1 



144 



Yi' 



Cg 



1 1, 



7^1 



18 



Yi- 



c.l-^-Y, 



25 25, 



:Yi' 



cs[l--Y, 



25 o 
144 ^ 



17, 



+^^ ' 16 - 24^^ + 72^^ + 384^^^^ 



29,, 2 1 



1 



1152 



Y1Y2. 



1 



331776 



+4 Y2+ce[-~-Y, + —Y,' + —Y2 - -— n Y, 



+C5 



33 

16 



96 12 

17 



144 



'-Y, + ^Y,^ 

4 144 1152 



^ 864 



1 



331776 



(i^2. 



(i^2, 



1728' 



tY 



3,2 



20736 



nis. 



C5 C6 ( 24 - -^2,1 + ^2^3,2 



55 _ J_ 1 

+^^' 4 432 ^-' + 20736 



Yi,: 



(4.4) 



The analytical form of w, 



1 — Symanzik ■ 



is very lengthy and its explicit expression is reported in Appendix D. 



Numerically these expansions read (we only report five significant digits, which are usually enough for simulation 
purposes, although the results shown in Appendix A for the finite integrals can be readily calculated with greater 
precision) 



E' 



■standard 



= 1- 



n-l 



9 



n-1 
^32" 



9 



0.0072699 (n - 1) + 0.0059930 [n - If 



9 



0.0028167 {n-l) + 0.0034299 {n - if + 0.0015673 (n - \f 



771O— Symanzik 

~ 12 



15 n-1 



g - 0.024449 {n -1) g^ ~ 0.0044905 (n - 1) + 0.0042082 [n ~ 1)' 



0.0014508 {n-l) + 0.0017541 [n - 1)^ + 0.0010241 (n - 1)^ 



771 1 — Symanzik 



15 n-1 



9 



(0.024449 + 1.60443 cg + 1.02179 cg) (n - 1) 



9 



12 4 

- (0.0044905 + 12.5286 c^ + 0.26629 cg + 4.7831 c^ + 0.44431 cg + 15.0482 cg cg 

+0.44752 cj + 0.36265 cg + 0.15246 cg) (n - 1) 



(0.0042082 + 0.68841 07 + 0.34420 cg + 0.050819 cg) (n - if 



9 



(0.0014508 + 0.11759 cg + 4.7668 Cg + 108.382 4, + 0.068884 cg + 5.8854 cg cg 

+ 186.142 c^ Cg + 1.7585 c^ + 111.566 cg c^ + 23.6332 c| + 0.11683 C7 
+9.2935 C5 C7 + 5.7226 cg C7 + 0.078923 cg + 7.1066 C5 cg + 4.5248 cg cg 
-0.031694 Cg + 2.0109 C5 cg + 1.2889 cg cg) (n - 1) 
+ (0.0017541 + 0.068786 C5 + 0.047639 cg + 0.15515 cy + 12.853 C5 cy 
+8.18909 Cg C7 + 0.069749 cg + 6.4265 C5 cg + 4.0946 cg cg 
-0.011464 eg + 0.67029 C5 cg + 0.42964 cg cg) (n - 1)^ 



+0.0010241 (n - l)'^ 



9 



(4.5) 



Notice that at order 0{g) the result is the same in all three cases. This is a consequence of the equipartition of the 
energy. The numerical value of the quartic coefficient in ij^tandard (jiffg;^-s roughly by 3% from the result in Ref. [|| . 
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This is due to the more accurate determination of the integral H^ obtained in the present paper (see Appendix A). 
This difference is too small to change at all any conclusion of Ref . Q . 

Now the energy-based effective coupling constant gs is defined in a non-perturbative way, 

9E = , (4.6) 

Wi 

where £''^^ is the Monte Carlo measured value of the internal energy at some value of the bare coupling g. The 
perturbative expansion of g^ in terms of this coupling g is obtained from the previously calculated expansion for E 

9E^9 + — .g' + — .g' + — / + ■ • • . 4.7 

Wi Wi Wi 

The lattice beta function (3^ in terms of the bare coupling g is |g^ 

/3^(.g) = ~a^ = -/?off' - /?i5' - P2g' - P39' -■■■ , (4.8) 

da 

and in terms of g^ it becomes 

2 o 3 02wf- Pi W1W2 - Po (wj - Wi W3) 4 

= ~Pa 9e - Pi 9e 2 ^ 9e 

/^s ^1 ~ 2 P2 wl W2 + f3i wi W2 + 2 /3o (2 wf — 3 wi W2 w^ + w\ 1^4) 



9%---- , (4.9) 



where the function 9{gE) is obtained by inverting Eq.(4.7). 

The integration of the beta function yields the dependence of the lattice spacing a on the coupling constant 

aA = {Po9r^'^^° exp (^-^^ (1 + 0(5)) , (4.10) 

where A is the integration constant, the lattice Lambda parameter. An analogous equation can be derived in the 



effective scheme with an integration constant A^. From Eq.(4.7) and Eq.(4.10) the ratio of these two constants can 
be exactly determined, 

Ab = Acxp(^^| . (4.11) 



V. CONCLUSIONS 

We have calculated the perturbative expansion of the internal energy for the tree-level and 1-loop improved 
Symanzik actions on the lattice for the 2D nonlinear cr~model with symmetry 0(n) up to fourth order in the coupling 



constant. Th ese r esults are shown in Eq.(4.5). The definitions that we have adopted for the internal energy E are 
shown in Eq.( 2.1C ). These expansions allow the definition of an effective coupling gE- We expect that the series which 
determine the lattice spacing in terms of the coupling constant are better behaved if expressed in powers of gs- This 
hope will be checked in a future publication where we plan to calculate the mass gap of the model through a Monte 
Carlo simulation. To this end, we have also given the analytic expression of the lattice beta function in terms of gE- 
The calculation is rather involved and it has been done separately by the two authors. Only the final results were 
compared. Also the numerical value of the integrals in Appendix A has been obtained independently and checked 
afterwards. Further checks were done: for example some subsets of diagrams must altogether yield an IR-finite result. 
These tests come out when we consider other definitions for the energy operator different from the one shown in 
Eq.( ^.lO ). For instance if we use the whole expression of the 1-loop Symanzik action as an energy operator, Eq.(^.5D 



then diagrams 9, 11, 12 and 13 of Figure 6 are multiplied by a different factor from the rest of diagrams. Then this 
subset, taken separately, must produce an IR-finite result. 
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VI. APPENDIX A 

List of finite integrals 

We will give the list of finite integrals that we have used to express the results for E as a manifestly finite quantity. 
The basic notation is introduced in section 3. Besides, the following definitions will be needed: 

2 

Ap,, = {p + q) -]? -q^ , 

\,q — ^P+q - lip - Eg , 

A'^ = TT — TT — TT 

A° , = □p+, - Dp - □, , 

A;^^^{pTq)l-pl-ql. (6.1) 

In the three-loop integrals we use the notation 



+" d2p /•+^ d^q [+'' d^fc /•+^ d^r ,,;,., , , 



-X (2t:YJ-, {2TifJ-, {2T:fJ-, (27r)''^"' 



and the definitions and numerical values of the integrals are 

^p^q ^k,r 

pp (p k^ r^ 



K= I D3 P,q ^k,r ^ 0.0958876 , (6.3) 



2 



J= I D3 ( Pi^) =0.136620, (6.4) 

-pA qZ ^A y,Z 

K^= I D3 ^ ^'"'''r, = 0.0673313 , (6.5) 



Hp n, Uk Ilr 


= 0.0673313 


[^pqkr) 


= 0.104551 , 


Up Uq Uk Ilr 



J"" = D3 ^ ^ = 0.104551 , (6.6) 



RS = I D3 ^ — '^^ — = 0.0572726 , (6.7) 

' {Hpf n, Uk Ur 



[^pqkr] P 

J'? = / D3 — ^ i = 0.0867807 , (6.8) 

(Upf Uq Uk Ur 

KS= I D3 "'" '''" = 0.0578002 , (6.9) 
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J'^ I Ds 't\'"Z - 0.0809553 , (6.10) 

lip iiq ll]^ li^ 



'^^ I ^^ §fB^r -'■'''''' ^ ^'-''^ 



52 = / As ^ = 0.184636 , (6.12) 



'^^ I ^^ ufu^t^lir -'■''''''' ^'-''^ 



S^ = I L»3 P"'"' ^ ^ = 0.148440 , (6.14) 

Ds ^t'^J''^' t''' - 0.261935 , (6.15) 

iip iig life llr 



S,^ I D, ^^^^^ ^"'"^ = -0.181963 , (6.16) 

Lip Liq Lif^ iij- 






S,o ^ I D, "^+^^\^ ^^-'- = 0.0881094 , (6.20) 

lip llq Ilk i-i-r 

Si. s / D, "";; ^'t,"J^'-' = 0.0649188 , (6.22) 

lip llq iLf^ iLj- 

5i3 ^ / D, ""t! n'^'n^'"' = 0-113020 , (6.23) 

lip llg 11^ 11^- 
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S,6 s / Da "';,' „"■„' „'■"' = 0.233959 , (6.26) 

lip ilg 11/^ II7- 

•^iT ^ / D, "p+^^p^^^^-^ ^ 0.0963352 , (6.27) 

iip iiq lifc iir 

?i8 = / £'3 \ = 0.184379 , 6.28) 

{Hpf Ug Uk n, 

f Af„Af^n 
5i9 = / D3 ^ — k^r_^_ ^ 0.120705 , (6.29) 

520 = / D3 ^ I ' = 0.335775 , 6.30) 

y (Hp)' n, Hfc n, 

/■ ( ^pqfcr ) ^P 

521 = / Dz —^ — = 0.213243 . (6.31) 

J {Hpf Ug Uk Ur 

The above listed integrals are not all independent. In Appendix B we give a few identities that these integrals 
satisfy. 

The measure for the four-loop integrals is 

, "^ d^p [+'' d^q /•+^ di^k [+"" d^r [+"" d^s [+'' dih 
Da = 



^ (27r)^ J^^ {2Trf J -^ {2Trf J -^ (27r)^ J ^^ {2T:f J -^ (27r)^ 

X {inf 6^{p + q + k + r) (27r)^ 5^{k + r + s + t) , (6.32) 

and the numerical values of the four-loop integrals are (the result of H3 differs in the third significant digit from the 
less accurate result given in ||]) 

H,= [da ^p^^ Afc,. Sp,^ ^ 0.0378134 , (6.33) 

J p2 qi ^2 ip. 'p. £2 

^^^ r^ Afc,r Sp,fc,- ^pgst ^ _o.0322778 , (6.34) 

J fp q^ k'^ r^ 'P t^ 

i?3 = [04 ^'^ J''" — "^ = -0.0128736 , (6.35) 

J p2 q2 f,2 p g2 £2 

jj^^ f j^^ ^pqkr ^krst Sp^t ^ 0.0411085 , (6.36) 

J p^ cp k-^ r^ s^ 1r 

H5= f Di ^p^g^fc.'- ^M ^ -0.0501531 , (6.37) 

J p2 qi j/2 ip. p, £2 
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^s ^s Y^s 
Hf = I D4, — — '"'"' '"^'^ = 0.0218345 

1 ' UpUgUkUrUsUt 



(6.38) 



m 



Da 



AS yS yS 
k.r pqkr pqst 

Hp n, Uk Ur Us Ut 



= -0.0181139 



(6.39) 



m 



Da 



A\ A^ A^ , 

p,k r.s q, — t 

n„ n„ Hfc n, n, n^ 



-0.0042338 



(6.40) 



y^S y5 y^S 

Hf = I D4 „ pgfc'-^fe';^* JPi^l = 0.0262036 



Up Ilq Uk Ur Us Ut 



(6.41) 



m 



Da 



AS ^s ^s 

p,q k,r s,t 



Up Ug Uk Ur Us Ut 

On the other hand the one-loop integrals are defined as 



= -0.0327709 



Y, 



-TT i2 



d'q fa 



Y ■ = 



+^ ^2 



-. i2nf \UJ ' "''' J_^ i2Trf{Ugy 

and their results are listed in Table 1. Notice that Ya =Yi. 



(2* > j) 



(6.42) 



(6.43) 



Table 1: One-loop integrals. 



n 


2.0435764382979844236 


Y2 


4.7830710733439886212 


Y3 


11.816615246907788250 


yi,2 


0.4729502261432961899 


^2,1 


30.077096804291341057 


Y3A 


558.65986413777280387 


Y3,2 


77.324121011413132160 


>4,1 


11817.841483609309517 


14,2 


1489.1480965521674895 


^4,3 


202.26364872706189510 


>5,2 


32200.496224041766111 


15,3 


4006.2729031961906982 


^6,3 


88276.902118545681915 



Some of the one-loop integrals were introduced in Ref. g. All li.j with i ^ j show up only in the results for the 
1-loop Symanzik action. 

In the perturbative expansion of the 1-loop Symanzik action there appear some vertices with a very high mass 
dimension. Once these vertices are inserted in the corresponding Feynman diagrams, all propagators may cancel 
leading to non-fractional integrands. In this case the following expression can be useful 






2m- 1 



m 



TO > 1 



For example. 



-TT i2 

d q 



. (27r) 



9' =4, 



J-^ (27r 



2 □, = 12 



+ TT 



(2^) 



jUg = 5 



(6.44) 



(6.45) 
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VII. APPENDIX B 

Identities among the integrals 

Some of the integrals Si in the above Appendix are actually related among themselves. We have not taken advantage 
of these relationships in the final expressions of section 4 because we are not sure to have discovered all of them. In 
this Appendix we show the identities that we have found out, 

l^Ss=K^-K^, ^S,, = K^-K^, l52i = J^-J^, (7.1) 



2 



2S +S I 5 - ^9 35^1^ , /65 5Y,,,\ Y,., jY,,,)^ 83 lir^g (^2,1) 

" ^° ^ 8 864 ^ V96 41472/ 576 165888 M 32 6912 1990656 



, (7.2) 



5*1 — 5*10 — 453 — 25*7 + -5*4 = 

67 , 41Y,' , 2 / 169 7^24 \ 7^2.1 (^2.1)' (^24)' ^39 ^2.1 (y2,i)' \ ,7 3^ 

1 I 1 no con/(/f I Aar\o cgocco oocGC/i/fC/f 1 I ic 1 'roc nncooc I ' v • / 



32 432 V 192 82944/ 4608 663552 286654464 \ 16 1728 995328 

We give the proof of the relationship Eq.( |7.2| ). Twice the numerator in Sn plus the numerator in 5*10 is equal to 

^p,q ^k,r (Hp+fc + Ilp+g + Hp+r) , (7.4) 



which, by using Eq.(3.5) under the integration, becomes 

Ap,,Afc,,(4np-E^^fc,) . (7.5) 

Therefore the following expression is true 

2Sn + S,o^^fD, ^M^_^,. (7.6) 

J ttg il/c iir 

The numerator in the integrand can be rewritten after some straightforward algebra under the integration 

Ap,, A,, = -i (k^y r' + \f (fcTry + 1 (f)' 1??^ , (7.7) 



and now a use of the second identity in Eq.(|3.6|) yields the final result shown in Eq.(7.2). 



VIII. APPENDIX C 
Methods for the numerical calculation of the integrals 

We have used three methods to calculate the various finite integrals listed in Appendix A. They are: i) an extrap- 
olation to infinite size of the results obtained at small sizes, ii) the Gauss integration and in) an extension of the 
coordinate space method p3| with the Symanzik propagators. Several integrals were evaluated by using more than 
one method for checking purposes. In this Appendix we will briefly describe methods i) and in). 

In the first method we calculated the integral / at small lattices L obtaining I{L) for several L. Then we extrapolated 
this set of results to infinite size. The extrapolating formula was 

, s , > &1 ^2 log L 

I{L) ^/(i^^) + -J- + A_^, (8.1) 

where the correct result for the integral is I{L — cxd). The exponent in the denominators is m = 1 for the one-loop 
integrals and m = 2 for all other integrals. The errors were determined by looking at the stability of the figures 
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against the increase of L. In fact such a stabihty indicates that further terms in the expansion Eq. (8.1) are irrelevant 
for the last stable digit. For one-loop integrals we obtained actually many significant digits by working up to lattice 
sizes as large as L — 10000. For the Si integrals the largest size was L = 60. 

For finite L the integral I{L) is actually a sum of terms. In this sum we have excluded the momenta which lead 
to vanishing propagators (jp = or lip = 0). In the limit L — > oo this procedure implies the exclusion of a region 
that has zero measure in the corresponding integration and therefore it has no consequences on the final extrapolated 
result I{L = cxd). We have checked this statement by repeating the calculation on several one loop integrals firstly 
i) by excluding only the zero mode j? — and secondly ii) by excluding the lines pf = or p^ = 0- Both methods 
yielded exactly the same final extrapolated number. 

The integrals Hi and Hf emerge in the calculation of the diagram 10 of Fig. |3[ They are four-loop integrals and 
as a consequence a direct application of the above method is rather slow and a poor precision is obtained. In order 
to simplify the evaluation of these four-loop integrals, one can take advantage of the topology of the diagram and 
rewrite them as an effective two-loop integral 

H. ^r^T. (fi r -^^"''"'f'l , m) 

and analogously for H^ . The sum ^j, over numerators TV/"' (p, q — p) contains one term in H^ and 64 terms in H^ . 
Nevertheless, in all cases the integration of Eq.(|8.2|) is much faster than a direct inte grat ion of Eqs.(|6.33)-(6.42). The 



largest lattice size used for the evaluation of these integrals with the method of Eq.( p.2D was L = 400. 

In some cases we have also applied the coordinate space method to check our results. We need to extend this 
method, introduced in p3| , p4| , to include the case of improved propagators. Here we will show with some detail the 
calculation of J . This requires the previous evaluation of the improved free propagator G{x) 



■" d^B e*''^ - 1 



In terms of G{x) the integral J^ reads 



-^^ - ^ E f 16 E (d^^d^Gi.))' 2 Y: {di {dt + d-) G{x)t + ^ E ((^^^ + d^) iK + d-,) GixA 



(8.4) 



where 9^ have been defined in Eq.(2.4). The procedure is analogous to what is done for the standard case, the only 



new features arise because of the presence of an improved Symanzik propagator. This propagator, Eq.(B.3), satisfies 

CG{x)^-5^{x), (8.5) 

where (5^(x) is the Kronecker function on a two-dimensional lattice and C is the improved Laplacian 

^-^E^X--i^E(^.^+^;)'- (8-6) 

If x^OEq.(|l|) provides 

E G(x + 2/2) = - E G{x - 2/2) + AG{x) + 16 ^ [G{x + /2) + G{x - /2) - 2G{x)] . (8.7) 

fi ^J. fj, 

Now, by using 

■ lognp = -pj- ( ^ sinp^ - - sin2p^ ) , (8.8) 



we obtain 



where 



dp^ ^ ' Hp V3 "^ 3 



- (G(.T + /2) - G{x -n))-7: {G{x + 2/2) - G{x - 2/2)) = x^ X'{x) , (8.9) 

3 6 
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X'^ix) ^ 



-" d^p 



logHp 



Eq. ( ^.9[ ) provides a recurrence relation for the propagator, 

G{x + 2/2) = Gix - 2/2) + 8 [G(a; + /2) - G(a; - /2)] - 6a;,,X^(x) 
which must be complemented with the relationship 



X^{x) 



E 



i-(ioc(.).EG-( 



X - 2/2) - 4 G{x - /2) - - G'(x + /2) 



(8.10) 



.11) 



.12) 



obtained from Eq.(8.7) and (|8.9|) a fter s umming over /i. 

From the recurrence Eqs.( pTll[ )-(8.12) and the symmetry relations G(a;i,a;2) = G{x2,xi), G{xi,X2) ~ G{—xi,X2) 
{xi and X2 are the two components of the coordinate site x) we can determine G{x) for any x in terms of the values 
of G{x) at a basic set of 10 sites. These sites are shown in Fig. ^ Below we give the values of G{x) on these points 



in terms of K^.j and Yj 



G(0,0 
G(1,0 

G(l,l 

G(2,0 

G(2,l 

G(2,2 

G(3,0 

G(3,l 

G(3,2 

G(3,3 



48 



Yi 



1 



1 

4 

-Y2,i - 



1152 
3 



1 1 

12^^^8 



27648'"^' 576 ^'^ 48 ^ 576' 

lA 1 "~ 1"^ 1 ~r y'? 1 — — j'^i 

331776 ' 1728 ' 144 ' 3 

13824 ■^' 32 ' 16 288 

^Y,,-^Y2, + ^Y,-'-^, 
1728 ' 128 ' 12 72 ' 



1 



lY^A - 



1 



917 
576 



:Y4 



7962624 ' 55296 
95551488 ■ 331776 



4,1 



27648 
1 



Y3,l 



—Y21 
64 ' 



119,^ 233341 

Yi H , 

48 55296 ' 



1 7 8Q 27 

y.5 1 + ^— n 1 Y3 1 + Y2 1 Yi 

'' 3072 ' 864 ^' 128 ' 4 



169681 
497664 



(8.13) 



A 



FIG. 7. The basic lattice sites for the evaluation of the propagator in the coordinate space method. 



We must know the values of G on the above set of sites very accurately but the above results for Yij are not precise 
enough (see Table 1). This goal is achieved by using the recurrence Eqs.(3.11)-(8.12) backwards. 



In some cases, we also need the calculation of the squared propagator G2\x) defined as 
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G2{X)= -^ ', ^ ^ ^ (8.14) 

-^ (2^)' (np)2 



The quantity TO has been defined in Eq.(3^). The laplacian Eq. (|8.6| ) apphed on G2 produces CG2{x) = —G{x), or 
alternatively 

Y^ G2{x + 211) = -Y^ G2{x - 2/1) + 4G2(a;) + 16 ^ \G2{x + Jl) + G2{x -'p)- 26*2(2;)] + 12 G{x) . (8.15) 

MM M 

Now, by using that 

we obtain the basic recurrence relation 

G2{x + 2/2) =G2{x- 2ii) + 8 [G2{x + fi)- G2{x - V)] + 6x^Xf (a;) , (8.17) 



where Xj (a;) = G(a;) + l/47r. Summing this equation over ^ and using Eq.(8.15), we get 



i ^G2(x-2m) + ^ ^G2(a; + /2)+4^G2(a;-/2)] 

M M ^1- J 



Xl{x) = ^^ [2G(:e)-10G2(x)-^ >^G2(x-2/2) + ;: > ^62(2; + /i) + 4 >^G2(a;-/2)) , (8.18) 

2^At^M \ 



which must be complemented with Eq.(B.17) to obtain G2{x) for any x from the same basic set of points shown in 
Fig. |. 

Further tricks necessary for the computation of the integrals are analogous to those already explained in Ref. [p3| . 

IX. APPENDIX D 

T71 n 1 — Symanzik 

Form or w^ 

The expression for w^ ymanzi ^^^ ^^ written as the sum of w^ ymanzi pj^^g ^ pure 1-loop Symanzik action 
contribution. This contribution contains terms proportional to (n — 1), (n — 1)^ and several powers of the coefficients 
Ci, i = 5, ... , 9. Let us parametrize these contributions in the following way 

+ o'*^ " ^> V"^5 95 + C6 ge + c^ qj + cg gs + cg gg + 05 cy 957 + C5 cg gas 

+C5 Cg (J5g + C6 C7 ge? + Cq Cs (?68 + Cg Cg (Jeg , 



1 — Symanzik — Symanzik 



(n-l) [C5 P5 + cg P55 + C5 P555 + Ce P6 + Cg P66 + Cg P666 

+c. 

+C5 C7 P57 + C5 C8 P58 + C5 Cg P59 + Ce C7 P67 + Cq Cs P68 + Ce Cg P69 ) • (9.1) 



Now the several coefficients qi are 

<l5 = x-Sis ^ -7S1 — S'2 — -5*3 , 

«7--i-i^i" + ig^i-^^i +6912^^ ^72^^ +144^^ ^^"10368^^ ^' ' 

.8 = 4- 1^12 + Si^i - T^n^ + 4Ln^ - ^y^ + in r^ - :ir,^ Y2 



4 96 1152 13824 144 288 20736 

'^^ = " 2 " 16 ^^" + '^'' ^^''^ + 3^^ 144^^' + 432^^' " 36^' + 48^^ ^' 1296^^' ^' ' 
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.5v = 8 - \y, + Ir,, - ^r, Y,, - ^ra. + ^n ^3. , 

«58 = 4 - -Fi + -Y2.I - — Fl ^2,1 - ^Y,^, + j^^^Y, ^3.2 , 

959 = -2 + -Fi + -^2,1 - 24^1 y^.i - 2Y6^3,2 + -^Y^ Y^.^ , 

967 = 4^1 - ^>1' - i>^2 + ^Yi Y2 , 
3 3 3b 

6 6 72 

969 - 2^1 - ^Fi^ - 1^2 + ^Fi ^2 , (9.2) 

and those pi are 

123 ^Ic Ic^lc^c^lc 1q^1q^q^q^93 185 , 215 3 

^^ = -128 + 6^^°-2^^«+4^^+^^+2^^-3^^+4^'^ + ^« + ^^+64^^-288^^ +2304^^ 

I 353 121 , 655 .y ^ ... ^2 7 ^2 1 . .3 
18432 ^'^ 6912 ^ ^^165888 ^ ^'^ 884736^ ^'^^ 1327104 ^^ ^'^-^ 1146617856^ ^'^^ 

5 145 175 2 5 25 11 

1^2 H Fi y2 Yi^ Y2 H 1^2 Y2 1 Yl ^2 F2 1 Y3 2 

48 1728 10368 10368 ' 124416 ' 6912 ^'^ 

77 605 2 1 5 

H Yl Fs 2 Yl Y3 2 + Y2 1Y32 Yl Y21Y32 

55296 ■ 1990656 ^'1990656 ' ' 23887872 ■ '^' 

5 25 1 5 25 

"248832^2 ^'^' + 2985984^^ ^' ^'^' + 20736^^^^ " 124416^^ ^*^' + 2985984^^' ^^^' ' 
557 425 349 25 7 2 13 5 
P"' - "1^ + ^^1 + T^^^'i - lii^i ^^'1 - ^^ (^2-i) + 1728^^'^ 1728^^ ^^'^ 

1 5 35 5 5 2 11 

2 9 1 y t 1 — 1 'i 9 ~h 2 1 2 S 9 "f' 29 1 2 ■^ 9 — ( 2 '^ 9 } — lA 9 

248832 ' "^^ 288 ' 5184 •^' 62208 ' ' 2985984 ^ ■^' ' 27648 ' 
5 1 15 

+ 31104^' ^^'' + 11943936 ^'-^ ^^'' + 124416^''' " 1492992^' ^''' ' 
_ 1670 _5_ 1_ 1 ^^ 

^'''" 27 +5184 "'' 31104 '''+2239488 ■' 

5 lo .s 7? lo Ir. 41^ 95 ^. 1889^^3 5^^ 

P6 = \- -S21 - J^ +J^ 5i9 + -58 H Yl Yl H Fi Y2 

^^ 8 6 2 '^ 4 ^ 32 128 13824 32 

29 55 n 5 n 25 2 1 5 25 2 

+ Yl Y2 Y2Y12 Y2' + Yl Y2' + Y3 Yl Y3 + Yi^Ya , 

192 1536 1728 20736 144 864 20736 

P66 = 6Y1 - ^ (Yi^ + Y2) + |yi Y2 - Ay,2 + ly3 _ ^nY^ , 

_ 4 

P666 = ^Ys , 

p.e = -15 + 12Yi - ^Yi^ + |y2,i - |yi Y.i - ^ (Y.i)' - ^Y2 + |yi Y2 + A Y2 Y.i 

-i^Y32 + ^YiY32 + ^Y2iY32-^Y2Y32 + ^Y43-^Yin3 

576^'' 3456 ' ''' 82944 "'' ''' 10368 ' ^'' 432^'' 5184^ *'' ' 

P556 = 44 + - Y3,i - — Y4,2 + ^Y^Y5,3 , 

2 1 

P566 = 4Y2.1 - -Y3,2 + :jtY4,3 , 

3 3d 

145 2y 1 n- 1^ I 11 y ry i^ ^^ y + ^9 ^ ^ 85 .^ 

"165888^^ ^^'^^ 82944^^'-^^ + 1990656^^ ^^^'^^ " 288^' + 384^^ ^' ~ 5184^^ ^' 
1 5 115 

2 9 1 2 9 A- 2 1 2 9 2 9 1 -\- 2 ■^ 9 —- 2 1 Y'i 9 "h 2 1 2 ■^ 9 

13824 ^'^ ' 165888 ' 4608 ^' 6144 ^ ^'^ 165888 ^' 
1 5 

+ 1990656 ^'^ ^'^ " 23887872 ^ ^'^ ^'^ ' 
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11 7 121 , 353 , 1 5 125 , 

^«-4-2^-+48^-576^^ +6912^^ -18^^ + 96^^ ^^-10368^^ ^- 

:Yi ^2,1 + TTTJTJ^^i^ ^24 - 777^2 + j^Yi Y2 - t^^^i^ Y2 , 



1728 ' 10368 ■ 12 16 432" 

_ 43 187 871 Vy + ^ ry ^2 ^Y ^ ^ YY 

P'' - " 16 + 1^^^ + 2304^^'^ 6^' ^'-^ + 3456 ^^'''> ' 576^^'^ + 1728^^ ^^'^ 

"248832^^-^ ^^'^ ~ 1728^^^' + 2592^^ ^'-^ ' 82944^^'^ ^^'^ + 27648^^'' 

"82944^^ ^^^^ + 11943936^''^ ^^'^ ' 

13 1 1 1 25 

P.S = -4 + -y, + -Y2, - -n r,! - ^^3, + ^n Fa. , 

P59 = -6 + 21-1 + Iy2^1 - ^Y, ^2,1 - ^y3,2 + ^n >^3,2 , 

-1 1 Oi^ 1 7 

Per = 9 + Yi - -Yi^ - —Y2 1 + — Yi y2 1 (y2 1)^ Y2 + — ^i ^^2 

- — y 1 J Q 9 -p ^9 1 J "^ 9 1 

576 ' 82944 ' ■ ' 

P68 = 4yi - ^^2 - ly2 + §yi y2 , 

3 3 36 

P69 = 6yi - 2yi2 - 2y2 + Vi y, . (9.3) 

Here a new set of integrals has appeared, Si, S2, etc., coming from diagrams 4, 5 of Fig. a. They are defined and 
calculated in Appendix A. This set of integrals is not completely independent and in Appendix B we give some 
relationships among them. 



.Yi-lYi^ 


48^^'^ 


576^^^^'^ + 


192'-^'^ 
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